Controllable generation of highly nonclassical states from nearly pure squeezed vacua 
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We present controllable generation of various kinds of highly nonclassical states of light, including 
the single photon state and superposition states of mesoscopically distinct components. The high 
nonclassicality of the generated states is measured by the negativity of the Wigner function, which 
is largest ever observed to our knowledge. Our scheme is based on photon subtraction from a nearly 
pure squeezed vacuum, generated from an optical parametric oscillator with a periodically-poled 
KTiOP04 crystal as a nonlinear medium. This is an important step to realize basic elements of 
universal quantum gates, and to serve as a highly nonclassical input probe for spectroscopy and the 
study of quantum memory. 

PACS numbers: 42.50.Dv, 03.65.Wj, 03.67.Mn 



I. INTRODUCTION 

Controllable generation of various kinds of nonclassi- 
cal states of light is crucial for the study of fundamental 
aspects of quantum mechanics and developing quantum 
information science. One of effective methods is a con- 
ditional preparation scheme by photon counting on one 
part of an entangled state produced in the parametric 
down-conversion (PDC) The main feature of the 

quantum states thus prepared is that their Wigner func- 
tions W{x^p), which is a quasi-probability distribution 
for non-commuting quadrature observables xp — px = i, 
should exhibit the negative values in certain phase-space 
regions. This is in sharp contrast to that a squeezed state, 
which are another representative of nonclassical states, 
have a non-negative definite Gaussian distribution of the 
Wigner function. 

Generation and observation of nonclassical states with 
negative values of W(x,p) have been reported recently. 
Those states can be categorized into two families. One 
family is the number states and their variants combined 
with coherent state components |^ ^ |^ , created in a 
non-colinear PDC configuration of the signal and idler 
photons. The other is the photon subtracted squeezed 
states, where a small fraction of a squeezed vacuum beam 
is beam-split and guided into a photon counter as trig- 
ger photons, and the remaining signal beam is output 
conditioned by the detection of the trigger photons [21 ■ 
These states must ideally be a particular superposition of 
odd number of photons, and can constitute superposition 
states of mesoscopically distinct components, sometimes 
referred to as the Schrodinger kittens Q ■ 

The negativity of the Wigner function can be easily 
lost with experimental imperfections. Among them, the 
imperfections of the homodyne detector are often cor- 
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rected to reconstruct a Wigner function for a generated 
state just in front of the homodyne detector. Direct mea- 
surement of the negativity of the Wigner function re- 
quires a generator with low loss and low noise as well 
as a nearly ideal homodyne detector. Negative values 
for the uncorrected Wigner functions were marked, e.g. 
W{0, 0) — —0.026 for the Schrodinger kittens in the 
photon subtraction scheme 0, and W(0,0) — —0.04 for 
the single photon state in the non-colinear PDC scheme. 
(The values correspond to a unit convention according to 
xp — px — i.) 

We present a scheme that can generate the photon sub- 
tracted squeezed states with the deepest negative dips of 
the Wigner functions ever observed to our knowledge, 
including the single photon state and the Schrodinger 
kittens. Realization of the single photon regime particu- 
larly requires a nearly pure squeezed input so as to collect 
appropriate trigger photons with high signal-to-noise ra- 
tio from a small fraction (10% at most) of the squeezed 
beam, which failed in the previous photon subtraction 
schemes. Our scheme enables one to vary the regime 
from the single photon to the Schrodinger kittens by sim- 
ply tuning the squeezing level which is directly controlled 
by a pump power for the squeezer, and hence to study 
continuous transition from the single photon regime to 
the mesoscopic superposition regime. 



II. EXPERIMENT 

The key is an optical parametric oscillator (OPO) 
with a periodically-poled KTiOP04 (PPKTP) crystal 
as a nonlinear medium, which can generate nearly 
pure squeezed vacuum states in a continuous-wave (cw) 
regime p, P. A schematic of our experimental setup 
is shown in Fig.^ which is similar to the previous 
work|^. A continuous- wave Ti:Sapphire laser (Coherent 
MBRllO) is used as a primary source of the fundamental 
beam at 860 nm, which is mainly used to generate sec- 
ond harmonic (430 nm) of about 200 mW by a frequency 
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FIG. 1: Schematic of experimental setup. OC:output 
coupler, TBS:tapping beam splitter, Triangle itriangle cav- 
ity, PBS:polarizing beam splitter, QWP:quarter wave plate, 
FP:Fabry-Perot cavity, HBS:50:50 beam-splitter. 



doubler (a bow-tie cavity with KNbOa), and is also used 
as a local oscillator (LO) for homodyne detection, and 
probe beams for various control purposes. The second 
harmonic beam is used to pump the OPO with a 10 mm 
long PPKTP crystal (Raicol Crystals) in an optical cav- 
ity (a bow-tie configuration with a round-trip length of 
about 520 mm). An output coupler (OC) of this squeezer 
cavity has a transmittance of 10.3%, and the intracav- 
ity loss is about 0.2~0.3%. The FWHM of the cavity is 
about 9.3 MHz. 

A small fraction (l'^10%) of the squeezed vacuum 
beam in path A is tapped at a beam splitter (TBS), 
guided into a commercial Si-APD (Perkin Elmer SPCM- 
AQR-16) through three optical filtering cavities in a row 
in path B, and is used as trigger signals for conditional 
photon subtraction. Figure[21 shows a photon counting 
spectrum of the trigger photons through the filtering cav- 
ities, which have 5~10 times wider bandwidths than that 
of the OPO, measured by detuning the OPO around 
the fundamental, degenerate, wavelength 860 nm. This 
consists of a single peak with a bandwidth of 8.6 MHz 
(FWHM), which can be easily matched with the observ- 
able mode by the homodyne detector. Other irrelevant, 
nondegenerate, modes from the OPO, peaking at every 
free spectral range of 573 MHz apart from the degener- 
ate frequency, are well suppressed sufficiently. The to- 
tal transmission ratio for the mode of interest is about 
30% just in front of the Si-APD. The trigger counting 
rate varies from less than 1 kcps up to 50 kcps, chang- 
ing with the pump power of OPO up to 160 mW and 
tapping ratios. These trigger counting rates are mostly 
much greater than the Si-APD dark counts (lOOcps). 
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FIG. 2: Photon counting spectrum of the trigger photons 
through the optical filtering cavities around the degenerate 
frequency as a function of cavity detuning: solid crosses for ex- 
periment, solid line for theory. The pump power was 160mW. 
All the OPO output is fed into the APD path for this exper- 
iment. 



The filter cavities are locked by a "sample-and-hold 
locking" technique which enables us to switch the system 
from a "locking phase" to a "measurement phase" period- 
ically. In the locking phase, the filter cavities are locked 
in a conventional manner (FM-sideband locking tech- 
nique). On the other hand in the measurement phase, 
the locking beam is cut off for photon counting and servo 
amplifiers hold the system in the same state as the one 
right before the locking beam is cut off. We attain this 
periodically to have the locking beam of the filters pass 
one optical chopper two times, before injected into filter 
cavities and after transmitted them before APD. Shut- 
tered timings are shifted each other in half period with 
respect to chopping frequency (500Hz). Uneven duty cy- 
cles of chopping disk prevent the locking beam from en- 
tering the Si-APD all the time. We also rebuild our servo 
amplifiers in order to be able to accept external timing 
signal and switch the two phases synchronized to the op- 
tical chopper's driver. 

The generated nonclassical light is combined with the 
LO at a 50:50 beam-splitter (HBS) and detected by a bal- 
anced homodyne detector with Si photodiodes (Hama- 
matsu S-3759, anti-reflection coated at 860 nm, 99.6% 
quantum efficiency). In order to improve the homodyne 
efficiency, the LO beam is spatially filtered by the mode 
cleaning cavity which yields the same spatial mode as 
the OPO output. The propagation loss (1^10%) mainly 
comes from the tapping itself, and the homodyne effi- 
ciency is 98 %. 

For every trigger signal from the Si-APD, a digital 
oscilloscope (LeCroy WaveRunner 6050A) samples ho- 
modyne signals over a period ~0.5/j,s around the time 
when the trigger signal is detected, and sends them to 
a PC one after another. Each segment of the homo- 
dyne signals are then time-integrated, being multiplied 
by a particular temporal mode function ^'o(i)i ^o provide 
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FIG. 3: Experimental Wigner functions (top panels) constructed from raw data without any correction of measurement 
imperfections in the case of 5% tapping ratio, (a) Single photon state generated by -0.7dB squeezed input, (b) and (c) 
Schrodinger kittens generated by -2.6dB and -3.7dB squeezed inputs, respectively. The insets in top panels are the contours of 
the Wigner functions. The middle panels are quadrature distributions obtained by homodyne detection. The bottom panels 
are photon number distributions obtained by the iterative maximum-likelihood estimation. 



the quadrature amplitude finally observed with respect 
to which the Wigner functions are constructed. About 
50,000 data points of the quadrature amphtude are col- 
lected, and the Wigner function is constructed usin g th e 
iterative maximum-likelihood estimation algorithm [lOj . 



III. RESULT AND ANALYSIS 

The temporal mode function 'l'o(i) should be chosen 
such that it defines the signal mode which shares the 
maximal entanglement with the trigger photon mode. 
In the Si-APD, a projection onto photon number states 
takes place in a time scale < 400 ps. The trigger outputs 
in this time scale are sent into the digital oscilloscope to 
sample homodyne signals. The time resolution in this 
sampling is about T — 1 ns. This is slightly longer than 
the Si-APD time scale, but much shorter than the tempo- 
ral correlation of the squeezed vacua whose bandwidth is 
typically 2B ^ 10 MHz. This T finally defines the trigger 
photon mode, which can simply be a rectangular tempo- 



ral mode function. For such a small i?T, a single mode 
description is valid [TH . In a good approximation, one 
can consider ^'o(i) in a formllJI (see also Appendix 

*oW = V^e-f°l*l, (1) 

assuming a trigger signal detected a,i t = 0, where 
Co = (7t + 7l)/2 determines the characteristic band- 
width Co/"' ^ 9.3 MHz of the cavity for the leakage rates 
7t ^ 57 MHz of the output coupler and 7^, ^ 1.2 MHz of 
the cavity loss. 

Figure|3| shows experimental Wigner functions (top 
panels), quadrature distributions over half a period (mid- 
dle panels), and photon number distributions (bottom 
panels) of the photon subtracted squeezed states. None 
of correction of measurement imperfections is made. 
From the left, the single photon state (Fig.Oi) by -0.7dB 
squeezed input, and Schrodinger kittens with two kinds 
of amplitudes (Fig.lJb by -2.6dB squeezed inputs andlS); 
by -3.7dB squeezed input). The tapping fraction of TBS 
in Fig. ^is set to 5%. The large negativity is obtained 
in a wide range of squeezing levels. 
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FIG. 4: The dependence of H^(0,0) on the squeezing level. 



Figure0]shows the values of the Wigner function at the 
phase-space origin W{0, 0) (solid circles) as a function of 
the ratio of the pump light amplitude to the threshold 
amplitude of the OPO, z = y/Ppump/Pth, which directly 
connected to the squeezing level in the cases of three 
tapping ratios 1%, 5%, and 10%. The lines are fittings by 
a theoretical model based on the two-mode state reduced 
to the signal mode \l/o(<) and the trigger photon mode 
00 (t) (Appendix El- 

The model includes two kinds of unwanted losses. One 
is in the homodyne channel modeled by a transmittance 
Th =78 %, taking imperfections of homodyne electronics 
into account. The other is a squeezing- level-dependent 
loss, Ts{z) — Tso — Kz'^, which also include a constant 
linear loss of passive optical elements. The latter is nec- 
essary to explain the degradation of the negativity of 
VF(0,0) as z increases. A loss induces a mixing of even- 
number of states into the photon subtracted squeezed 
state which ideally contains only odd- number states. 
Since the values of W{0, 0) for even-number of states are 
positive, the negativity of W{0, 0) is easily destroyed by 
a loss. The z-dependence of W{0, 0) cannot be explained 
by any model with constant linear loss. The optimal fit- 
ting is attained for Tso = 0-95 and k — 0.93. 

One of main causes for z-dependent loss might be phase 
fluctuations of the LO in the homodyne detector and of 
the locking beam at several cavities. Such fluctuations 
cause effectively phase noises which becomes larger for 
larger radius in the phase space. As z increases, the 
Wigner function spreads over the phase space, and suffers 
from such fluctuations. For realizing larger Schrodinger 
kittens, it will also be important to suppress phase fluc- 



tuations. 

In the single photon regime, i.e., small z, the trigger 
count rate becomes small. A primary factor to determine 
the behavior of W{Q, 0) is then the noise count rate v 
in the trigger channel, which include the dark counts of 
the Si-APD and background photons. It is necessary to 
reduce this rate to generate the single photon state with a 
deep negative Wigner function. In this regard, the OPO 
with PPKTP works very well to reduce noisy photons. 



IV. CONCLUDING REMARKS 

We have presented the photon subtracted squeezed 
states with large negative dips in their Wigner functions. 
Those highly nonclassical states cover the single photon 
regime and the Schrodinger kitten regime. Their contin- 
uous transition can be controlled by the squeezing level 
of the input source. The nearly pure input source from 
the OPO with PPKTP enabled us to study the degrada- 
tion of the negative Wigner functions in detail. It was 
found that there was the squeezing-level-dependent loss, 
whose origins remains open, but we suspect at present 
that it stems from phase fluctuations of the LO and the 
locking beam. 

Our scheme will particularly be useful to make the out- 
put states coupled to atomic media for quantum memory, 
and also to apply them to spectroscopic applications. In 
addition, our scheme can directly be combined with well- 
matured technology of Gaussian operations on quadra- 
ture observables in a continuous- wave scheme 

Mm for 

pursuing further challenges. For example, the attained 
level of the nonclassicality, i.e. the large negativity, will 
allow one to teleport the nonclassicality of non-Gaussian 
input states, which remains an unreached target in spite 
of its importance in quantum information science, by us- 
ing the presently available squeezer with -7dB squeez- 
ingiia. 

Another important chall eng e is to combine photon- 
number resolving detectors [To . llTl IT^ with our scheme 
for breeding the Schrodinger kittens into cats, and for 
developing a universal quantum gate circuit both for dis- 
crete and continuous variables [l3 . l20l l2ll| . 
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APPENDIX A: MODE FUNCTIONS 

To describe the trigger photon mode locaHzed in the 
time domain, an appropriate basis is given by 



1 , 27rfci, T T 

— exp(— i— -— ) <t< — , ,.,s 

T ' 2 - - 2 (Al) 

otherwise, 



In the present case where BT <C 1, the lowest mode 
is excited in the Si-APD with the weight ^ BT, and 
the weights of the other higher modes are negUgible • 
The trigger mode function is thus (f)a{t). The total quan- 
tum efficiency for this mode is given by 77 = rjorjpBT 
where rjQ is the intrinsic quantum efficiency and 77_f is the 
transmittance through the optical filters. Taking a noise 
count rate v further into account, detection of the trig- 
ger photons is described by the on/off detector POVM 
as specified in ref . ^ . 

The signal mode ^o(0 should be appropriately deter- 
mined by the statistics of the conditional output state pA 
in path A given by 



PA 



Tr ABipAsti^niV,'^)] 



(A2) 



where the state pab corresponds to the split beams in 
paths A and B just after the TBS (Fig. P). We consider 
a mode expansion of the field operator to describe the 
state PA 



K-l 



A{t) ^ 



fc=0 



where 



(A3) 



(A4) 



The mode functions {4'fc(i)} are determined by minimiz- 
ing the square error 



Jk = Tr 



.00 K-l 

/ dt A{t) - Ak^k{t) 

fc=0 



K-l 



dtTT[A\t)A{t)]- J2 Tr[4ifc].(A5) 



fc=0 



Since the first term is constant, we are to maximize the 
second term under the constraint of Eq. (|A4I) for each k. 
It is known that the optimal solution for this is given by 
the eigenfunctions of the integral equation 



Xfc*fe(i) = / dt'hit,t')^kit'), 



(A6) 



where 



h{t,t')=TrA[PAAHt)A{t% (AT) 
assuming that the trigger signal is detected at t = [l^, 

m 



In a good approximation, one can consider only the 
lowest order mode ^o(^) in a form 



*o(0 



<o\t\ 



(A8) 



The Wigner function is constructed with respect to the 
field component defined by 



A = Ao 



dtA{t)'^oit)- 



(A9) 



APPENDIX B: FORMULA OF WIGNER 
FUNCTION 

The Wigner function can be obtained by evaluating the 
covariance matrices in terms of the signal mode 5*0 (i) and 
the trigger photon mode <poit) as prescribed in ref. |12|. 
exploiting the on/off detector model in ref. [Tl| : 

W{x,p) = (Bl) 



l-AA(r;,z.) 



where 
R{x,p;r],iy) 



Mill, v) 



X exp|^ 

with 
and 



1 - T/i + T/iCT2(z) 1 - T,i + r,i(T2(_2) 



(B2) 



26-" 

V{2+[6(z)-l]r;}{2+[6(-z)-TM' 

(B3) 



(7^(z) = a{z) — 



(?(z)r\ 



where 
a(z) — 1 — T + T 



6(z) = T -f (1 - r) l-r,(z)7Tr 



2+[6(z)-l]r?' 

2Tsiz)^T z{z + 3) 
Co {z + l){z + 2y 



z + ir 



c{z) = v/(1-t)- 



- 2ts{z)^tVt z 

^ (z + 1)(z + 2)' 



(B4) 



(B5) 



(B6) 



(B7) 



Here r is the transmittance of TBS in Fig. ^ t/,. is the 
effective transmittance of the homodyne channel, and 
Ts{z) = Tso — Kz'^ is a phenomenological model for a 
squeezing-level-dependent loss, which is necessary to ex- 
plain the degradation of the negativity of W{Q, 0) as z 
increases. 



6 



[1] B. Yurke and D. Stoler, Phys. Rev. A 36, 1955 (1987). 
[2] M. Dakna, T. Anhut, T. Opatrny, L. Knoll, and D.- 

G. Welsch, Phys. Rev. A 55, 3184 (1997). 
[3] A. I. Lvovsky, et al., Phys. Rev. Lett. 87, 050402 (2001); 

A. Zavatta, S. Viciani, and M. Bellini, Phys. Rev. A 70, 

053821 (2004). 
[4] A. Ourjoumtsev, et. al., quant-ph0608230 
[5] A. I. Lvovsky, and S. A. Babichev, Phys. Rev. A 66, 

011801 (2002); A. Zavatta, S. Viciani, and M. Belhni, 

Science 306, 660-662 (2004). 
[6] A. Ourjoumtsev, et al.. Science 312, 83 (2006). 
[7] J. S. Neergaard-Nielsen, et.al., Phys. Rev. Lett. 97, 

083604 (2006). 
[8] T. Aoki et al.. Opt. Exp. 14, 6930 (2006). 
[9] S. Suzuki et al., Appl. Phys. Lett., 89, 061116 (2006). 
[10] A. I. Lvovsky, J. Opt. B: Quantum Semiclass. Opt. 6, 

S556 (2004). 

[11] M. Sasaki, and S. Suzuki, Phys. Rev. A 73, 043807 



(2006). 

[12] K. M0lmer, Phys. Rev. A 73, 063804 (2006). 
[13] H. Yonezawa et al.. Nature 431, 430 (2004). 
[14] N. Takei et al., Phys. Rev. Lett. 94, 220502 (2005). 
[15] S. L. Braunstein and H. J. Kimble, Phys. Rev. Lett. 80, 
869 (1998). 

[16] E. Waks, et al., Phys. Rev. Lett. 92, 113602 (2004) 
[17] D. Rosenberg, et al., Phys. Rev. A 71, 061803, (2005). 
[18] M. Fujiwara, and M. Sasaki, Opt. Lett., 31, 691 (2006). 
[19] D. Gottesman, A. Kitaev, and J. Preskill, Phys. Rev. 

A64, 012310 (2001). 
[20] S. D. Bartlett and B. C. Sanders, Phys. Rev. A65, 042304 

(2002). 

[21] M. Sasaki, et al.. Quantum Communication, Measure- 
ment and Computing, 44, (eds. S. M. Barnett et al., AIP, 
New York, 2004) ; quant-ph /0601058 



